Kinetics of ballistic annihilation and branching 
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We consider a one-dimensional model consisting of an assembly of two-velocity particles moving 
freely between collisions. When two particles meet, they instantaneously annihilate each other and 
disappear from the system. Moreover each moving particle can spontaneously generate an offspring 
having the same velocity as its mother with probability 1 — q. This model is solved analytically in 
mean- field approximation and studied by numerical simulations. It is found that for q = 1/2 the 
system exhibits a dynamical phase transition. For q < 1/2, the slow dynamics of the system is 
governed by the coarsening of clusters of particles having the same velocities, while for q > 1/2 the 
system relaxes rapidly towards its stationary state characterized by a distribution of small cluster 
sizes. 

PACS numbers: 82-20. Mj, 05.20Dd 



I. INTRODUCTION 

Ballistically-controlled reactions provide simple exam- 
ples of non-equilibrium systems with complex kinetics 
and have recently attracted a lot of interest They 
consist of an assembly of particles moving freely between 
collisions with given velocities. When two particles meet, 
they instantaneously annihilate each other and disappear 
from the system. 

Depending on the initial velocity distribution, two 
classes of asymptotic states have been observed in one di- 
mensional systems. I n g eneral, for continuous initial ve- 
locity distribution PJlfj|, as well as for some special case 
of discrete velocity distribution (symmetric two-velocity 
distribution ElUjq], or symmetric trimodal velocity dis- 
tribution with a sufficiently small fraction of immobile 
particles the steady-state turns out to be empty 

and it is approached algebraically in time. The dynam- 
ical exponent characterizing the time decay depends on 
the initial velocity distribution and it is still not com- 
pletely clear how to characterize the universality classes 
for this problem |^o| . On the contrary, for some dis- 
crete velocity distribution, the stationary state may not 
be empty, but may contain particles moving all with the 
same velocity (for example non-symmetric bimodal veloc- 
ity distribution [QJ^] or a trimodal velocity distribution 
with more than 25% of particles initially at rest 
This non-interacting state is generally approached with 
an exponentially fast decay. 

A richer behavior can be expected in a system with, 
in opposition to the ballistic annihilation case, an inter- 
acting steady-state. This can be achieved by constantly 
bringing new particles in the system by some suitable 
mechanism. A possibility is to allow branching processes: 



ballistically moving particles can spontaneously generate, 
with a given branching rate, some offsprings. Accord- 
ingly, one speaks of ballistic branching-annihilation. 

The problem of branching-annihilation has been re- 
cently studied in the framework of a diffusive dynam- 
ics [[ll],[l2] . The simplest example of such a system would 
be one with a single species of particle A, undergoing dif- 
fusive behavior, single-particle annihilation A — > 0, and 
branching A — > 2 A. There is always a trivial absorbing 
state, with no particles. For sufficiently low branching 
rate, this is the only stationary state, but for larger val- 
ues of this rate, another non-trivial 'active' stationary 
state appears. This stationary state phase transition be- 
longs to the directed percolation universality class [fl3|| . 
A slightly more complicated class of model are reaction- 
diffusion systems with the underlying reaction processes 
2A — ► and A — > (m + I) A, with m even. It turns out 
that for these models the critical exponents are not the 
ones of directed percolation but belong to a new univer- 
sality class characterized by branching and anni- 
hilating walks with an even number of offsprings. The 
constraint of local 'parity' conservation is the reason for 
the existence of this new universality class. 

Our aim here is to study the problem of ballistic 
branching-annihilation (BBA) in one dimension for which 
interesting new properties can be foreseen. The paper is 
organized as follows. In section |J, the BBA model is 
defined. The exact dynamical equations of motion are 
derived for the one dimensional case. In section III, 
the dynamics of the model is studied within a mean- 
field like approximation. In particular, the phase dia- 
gram of the steady-state is established in terms of the 
different parameters of our model. In this approximation, 
the steady-state is always approached exponentially fast. 
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Section IV is devoted to numerical simulations of the one 



dimensional model. It is shown that fluctuations plays 
a crucial role. Indeed, as in the mean-field approxima- 
tion, a phase transition occurs when the probability that 
the offspring takes the velocity of its mother is q = 1/2; 
however, for q < 1/2 the dynamics is be governed by 
the coarsening of clusters of particles having the same 
velocity, and the system approaches a completely filled 
stationary state with a power law decay. For q > 1/2, 
there is no coarsening and the system relaxes rapidly to- 
wards a non-filled stationary state. Finally, the results 
are discussed in section |v|. 



II. THE MODEL 

We shall first define precisely the BBA model stud- 
ied and secondly derive the corresponding equations of 
motion. 



2. A particle with coordinates (x, c) produces 
with probability pqdt a pair of particles with 
coordinates (x — a, — c ) and (x, c). 

3. A particle with coordinates (x, — c) produces 
with probability p(l — q)dt a pair of particles 
with coordinates (x, — c ) and (x + a, — c). 

4. A particle with coordinates (x, — c) produces 
with probability pqdt a pair of particles with 
coordinates (x, — c ) and (x + a, c). 

(the particular choice of the position of the newly cre- 
ated particle has been made in order that, independently 
of its velocity, a child cannot collide with its mother at 
birth.) Thus the parameter < p < oo characterizes the 
overall branching rate, while the parameter < q < 1 
characterizes the probability that the offspring has a ve- 
locity opposed to the one of its mother. The particular 
case p — corresponds to the pure ballistic annihilation 
problem previously studied |l|,||,||,|| . 



A. Definition of the model 

We consider a one-dimensional system composed of 
particles of size a initially uniformly randomly dis- 
tributed in space. Moreover, at t = 0, the velocities 
of the particles are random independent variables dis- 
tributed with the symmetric bimodal distribution: 



P{v) = -[S(v - c) + 8(v + c)] 
The dynamics consists of two mechanisms: 



(1) 



• The ballistic annihilation: Two particles making 
contact (with opposite velocities) disappear instan- 
taneously. 

• The branching: during the time interval [t, t + dt], 
the following branching processes take place: 

1. A particle with coordinates (position and ve- 
locity) (x, c) produces with probability p(l — 
q)dt a pair of particles with coordinates (x — 

a, c ) and (x, c). 
I 



B. Exact equations of motion 

We can now derive the equations of motion describing 
the dynamics of the system. In the particular case p = 0, 
a kinetic equation for the two-particle conditional distri- 
bution of nearest neighbors was derived as a rigorous con- 
sequence of the dynamics of ballistic annihilation |^,^|. 
This equation completely described the evolution of the 
system when initially higher order conditional distribu- 
tions factorized into products of two-particle ones. It 
was then possible to extract exactly and analytically the 
long time behavior of the particle density for several ve- 
locity distributions. Unfortunately, this property is no 
longer valid in the case with branching. Having not been 
able to find an observable in which one is able to repro- 
duce this exact closure, one has to face the usual problem 
of dealing with a complete hierarchy of coupled equa- 
tions |Q. It seems thus hopeless to find an exact ana- 
lytical solution to these equations. Accordingly we shall 
only write the equation for the one-particle density dis- 
tribution px(x,v;t). In section III, this equation will be 



solved using a mean-field approximation. 



A careful bookkeeping of the possible dynamical processes leads to the following equations: 
(d t + cd x )pi(x 7 c;t) = -2cp 2 {x 7 c;x + a, -c;t) 

+ P1 



px{x - a, -c;t) - y~] / dyp 2 (x-cr, -c;x + y,v,t) 

v=±c Jo 

P(l-g) Pi(x + a,c;t) - ^ / dy p 2 (x - y,v;x + a,c;t) 



V — ±C 



(2) 



and 



(d t - cd x )pi(x, -c; t) = -2cp 2 (x, c;x + a, -c;t) 



pq 



V — ±C 



pi(x + a,-c;t) - / dy p 2 (x - y,v;x + <r,c;t) 
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+ p(l - q) pi[x - a,-c;t) - 



where P2(xi,v\; X2,V2',t) is the joint two-particle den- 
sity to find a particle in the state (xi, v±) simultaneously 
with another in the state (X2, v^j at time t. 

The right-hand side of equation (0) can be interpreted 
in the following way: the first term describes the anni- 
hilation of a particle (x, c) with a particle of opposite 
velocity. It is given by the product of the density of 
a collision configuration [ap2{x, c; x + a, — c, t)] with the 
frequency of such an encounter (2c/ a). The second term 
describes the branching of a particle of velocity — c, at 
position x — a, giving birth to a particle of velocity +c 
at position x. This is only possible if no other particles 
are present in the interval [x, x + a] (otherwise there will 
be an overlap between two particles) and it happens with 
the rate pq. Finally, the third term describes the creation 
with rate p(l — q) of a particle whose mother have the 
same velocity. The same restriction as in the previous 
case applies. 

One can in principle write the equation of motion for 
P2 along the same lines. However, we shall not give here 
this cumbersome equation, as we are not going to use it. 

For simplicity, we shall only consider spatially homo- 
geneous system. We can thus write p\(x,v\t) — p\(v,t) 
and p2{xi,vi]X2,v 2 ]t) = p 2 {x\ — x%, V\\ 0, V2]t). Intro- 
ducing then the observable 'J'(i) = pi(c, t) — Pi{— c, t), 
one easily shows that it is an exactly conserved quantity 
when q — 1/2. This feature reflects the particular choice 
of rule, which are precisely symmetric when q = 1/2. 
As a consequence, one expects our model to exhibit a 
particular behavior at this point. 



III. MEAN-FIELD ANALYSIS 

A first attempt to obtain information about our model 
is to apply a mean-field approximation on equations (j^) 
and (H). One then assumes the following factorization: 

P2(xi,v 1 ;x2,v 2 ;t) = pi(x 1 ,vi;t)pi(x2,v 2 ;t) 

= Pi(v 1 ,t)p 1 (v 2 ,t), (4) 

(the last equality holds for a spatially homogeneous sys- 
tem) . 

It is then suitable to introduce in addition to the vari- 
able ^ the second variable: 



= pi(c,t) + Pl (-c,t) 
Equations (|j) and (||) lead to 



dt 



P(1-CT$)$-C($^-^), 



(5) 



(0) 



V — ±c 



dy p 2 {x - cr, -c; x + y,v; t) 



(3) 



d$> 

— =p(l-2q)(l-a^. 
The formal solution of this last equation is 

* (i) = *(0) exp \p(l -2q)(t-a [ dr $(r) 



(7) 



(8) 



As before, one sees that the value q = 1/2 plays a special 
role. Indeed, two regimes have to be distinguished: 

1. For < q < 1/2: the exponential term diverges 
unless (1 — (7$) — > as t — > oo. Thus a possible 
stationary solution is 



\Vl A — — 

* S — o j 

a 1 



(9) 



In the particular case q = 0, the time dependent 
solution can be obtained explicitly as shown in Ap- 
pendix. For t — > oo, one recovers the above sta- 
tionary solution. 

2. For 1/2 < q < 1: in this possible stationary 

solution is 



0, $ s 



-(l + c/pa)-\ 



(10) 



Is is straightforward to verify that the above stationary 
solutions are stable and are approached exponentially in 
time. 

Moreover, when q = 1/2 the complete time dependent 
solution can be obtained. From equation (0),one indeed 
finds ^(t) = const = and thus equation (g), becomes 



~dt 



p$ - (c + pcr)$ 2 + c<5q, 



(11) 



whose solution reads 

p 7 A cosh(At) + A(c + pa) smh(At) 



$(*) 



2(c + pa) A cosh(At) + 7(0 + pa) sinh(At) ' 

(12) 



with A = p 2 /^ + c(c+pa)^l and 7 = <&(0)-p/[2(c+pa)]. 
The stationary state is then given by ^ s = ^0 and 



$ s (g = 1/2) = (p + yp 2 +4c(c + pct)*§) j{c + pa). 

(13) 



and 



Here again, one sees from equation (|12|) that the steady 
state is approached in an exponential way. As already 
noted, ty(t) is an exactly conserved quantity for q = 1/2. 

The mean-field stationary phase diagram is shown in 
Fig. |l|. The stationary value $ s is plotted against q for 
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a fixed value of p ^ 0. The interesting feature is the 
presence of a gap A(p) for q = 1/2 given by 

A(p) = K 1_ TT^)' (14) 

A(p) decreases as p increases. When q < 1/2, $ s = 1/a 
for all values of p (completely filled state), while for 
q > 1/2, <f> s increases monotonically with p. The depen- 
dence is linear for small p, but <I> S — > l/cr when p — > oo. 

cr$ s 
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1 i 3&- 



I 1 
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FIG. 1. Mean-field phase diagram: the stationary value of 
the averaged density <E> S is plotted against q for a fixed value 
of p. 



3. Branching: for each remaining particle, one draws 
two random numbers, r p and r q , uniformly dis- 
tributed in the interval [0,1]. One offspring par- 
ticle is added to the left (right) nearest neighbor of 
a particle with velocity d (0) if the site is empty 
and if r p is the less than a given value p. Hence, 
p is the probability of branching. This new parti- 
cle takes the velocity of its mother with a certain 
probability 1 — q, i.e. if r q > q (and the other veloc- 
ity otherwise). If two particles are created on the 
same site (thus born from two different mothers), 
they annihilate instantaneously. 

For each of the above different steps, the sites were up- 
dated simultaneously. The simulations were run on a 
Connection Machine CM-200 and the data averaged over 
10 independent realizations. The mean initial density for 
all the simulations was 0.5, with, in average, the same 
densities of both kinds of particles. We have also shown 
that our results obtained for lattice of 2 17 sites were free 
of finite size effects. 

Note that when a particle branches, it can create at 
most one particle during a time step r. As a consequence, 
this limit the value of p that can be explored through the 
simulations. Indeed, the branching rate p is related to p 
via 

p = pr. (15) 



IV. NUMERICAL SIMULATIONS 

In view of the situation when p = one can 

anticipate that the fluctuations will also play an impor- 
tant role in the case with branching. One way to deal 
with the complete problem, including fluctuations, is to 
perform numerical simulations. 

The simulations were performed for a one-dimensional 
periodic lattice with typically 2 17 sites. The velocity of 
each particle was drawn from a symmetric bimodal dis- 
tribution. However, on computational grounds, the par- 
ticle velocities were chosen to be (0, c') (with c' > 0). 
The results for our model defined in section |l| can be 
recovered by performing a simple Galilean transforma- 
tion and putting c = c'/2. The particle size a is the 
lattice spacing, and the discretized time step is given by 
t = a jc! . 

The algorithm used to simulate the dynamics is the 
following. During one time step r, the three following 
processes occur sequentially: 

1. Ballistic motion: independently of the occupation 
state of the sites, the particles with velocity c' move 
one site to the right. 

2. Annihilation: two particles located on the same site 
disappear. 



Thus using the definition of r and c', one finds 



p being a probability, the adimcnsional branching rate 
pa J c can only take values between and 2. 

We can now discuss the numerical data obtained us- 
ing the above algorithm. Two kinds of quantities have 
been investigated: first, the time dependent density with 
particular emphasis on the stationary states and the way 
these stationary states are approached; second, a more 
microscopic quantity, namely the time dependent cluster 
size distribution P(l, t) in the system and some of its mo- 
ments. These quantities are well suited to describe the 
coarsening process present in the system. 

As in the mean-field approach and as expected from 
the last remark of section [n], the value q = 1/2 turns out 
to play a particular role and three regimes have to be 
distinguished. 

1. For < q < 1/2: The time evolution of the parti- 
cle density §(t) is shown in Fig. || for several val- 
ues of p. Clearly, the system reaches a stationary 
state $3 = l/u in agreement with the mean-field 
prediction. However, as shown in Fig. 0, the sta- 
tionary state is approached as $ s — $(t) ~ t~ x / 2 . 
This power law establishes after a crossover time 
roughly proportional to 1/p. 
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FIG. 2. Time evolution of the particle density a&(t) as a 
function of time t for q — 0.1 and several values of p. 

2. For 1/2 < q < 1: The time evolution of the particle 
density <&(t) is shown in Fig. || for several values of 
p. As depicted on Fig. ||, the stationary value of 
the density depends both on p and q. For p < 0.1, 
it is well fitted by 
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FIG. 3. <7$ s — a<b(t) versus t in a double logarithmic scale, 
for q = 0.1 and several values of p. For comparison, the full 
line represents t -1 ^ 2 . This decay establishes after a crossover 
time which behaves as r/p. 



4 = 0.9 



10" 



®s(p,q) «pexp(0.55/g), 



(17) 



Moreover, for p large enough $ s is not increasing 
monotonically as a function of p, but <I> S exhibits 
a maximum and then slightly decreases as p in- 
creases. As shown in Fig. |], the stationary state 
is approached in an exponential way according to 
$ s — $(t) ~ exp(— Apt), where A may depend on 

q- 

The limit case q = 1/2 is more difficult to inves- 
tigate due to the slow decay towards the station- 
ary state. In fact for p > 0.3, there are evidences 
that the stationary state is completely filled, i.e. 
$ s = 1/ct. For smaller values of p the simulations 
do not allow us to draw any conclusions, as shown 
in Fig. 0. 

Nevertheless, for q < 1/2, one has <I> S = l/cr, for all 
values of p while for q > 1/2, equation ([17]) shows 
that, at least for small p, <I> S ^ I /a. Thus for small 
p, $ s has a jump at q = 1/2, and we believe that 
such a jump will be present for all finite values of 
P- 
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FIG. 4. Time evolution of the particle density a&(t) as a 
function of time t for q — 0.9 and several values of p. The 
stationary state is reached after a time of order Wr/p. 
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FIG. 5. The stationary values of the averaged density er$ s 
is plotted against p, for several values of q > 1/2, obtained by 
numerical simulations. 
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FIG. 6. Semi-logarithmic plot of a§(t) — <r<& s versus t for 
q = 0.9 and p = 0.01. The exponential approach towards the 
steady state establishes for t/r — 250. 



10" 



4 = 0.5 

OOCBHIB) OOOSSXE) O 



o° ^ 







10- 



op = 1.0 
+ p = 03 
□ p = 0.1 
< p = 0.03 
A p = 0.01 



10" 



10' 



10" 



10' 



10' 



//x 



FIG. 7. Time evolution of the particle density a&(t) as a 
function of time t for q — 0.5 and several values of p. For 
small values of p (less than 0.3), we are unable to extract the 
steady-state density, for CPU reasons. 

We can now consider the properties of the clusters 
present in the system at a given time. The qualitative 
situation is well illustrated by the two snapshots in Fig. 

They represent the time evolution of a 512-site system 
during 1024 iterations. Moreover, a change of reference 
frame has been performed such that the particle veloc- 
ities appears to be ±c. Depending on q, one observes 
totally different pictures. In the case p = 0.7, q = 0.1 
(Fig. ||a), large clusters (of similar particles) are present. 
They are separated by two types of interfaces: vertical 
ones (which are stable) and rough ones. The dynamics of 
the system is totally governed by the random walks of the 
rough interfaces. During the time evolution, one rough 
interface may collide with a stable interface leading to 
the coalescence of two clusters into a large one. In the 
case p = 0.7, q — 0.9 (Fig. ||b), the sizes of the clusters 
are rather small and there is no stable interfaces. The 
dynamics is of a different type. 

A more quantitative description is given by the inves- 
tigation of the time dependent cluster size distribution 
P(£,t). In the domain < q < 1/2, where coarsening 
is observed, one expects |l5| that P(£, t) will obey to a 
scaling form: 



P{l,t) ~£- a n(£t-' 3 ). 



(18) 



G 





FIG. 8. Time evolution (vertical axes) of the configurations for a chain of 512 sites (the initial density is approximately one 
half) and for 1024 time iterations. The white pixels indicate sites without particle, the grey ones, sites with a particle towards 
the right and the black ones, sites with a particle moving towards the left. Fig. a is for p = 0.7, q = 0.9 while Fig. b is for 
p = 0.7, q = 0.1. 



In Fig. |9[ we plot the scaling function obtained by the 
collapse of the data for p — 0.7,q — 0.1, with a = 1 
and (3 = 0.5. Although the plot is very noisy, one still 
notes that the scaling function n(z) has a very particular 
shape, with a sharp maximum at z = z max . The value of 
Zmax increases slowly with q, going from 0.4 for q = 0.1 
to 1.2 for q = 0.4. 

A better way to extract the exponents a and f3 is to 
consider the n-th order moments of the distribution de- 
fined as: 



are of the form a n — —0.26 + 0.61n which is not com- 
patible with the relation ( po| ) . This probably shows that 
the simulations have not yet reached the true asymptotic 
regime. Moreover, as shown on Fig. [fl| P{£, t) is of the 
form: 



q = 0.1 p = 0.7 



10" 



(n 



f™d££ n P(£,t) 
Cd£P(£,t) 



(19) 



which according to the scaling form given by equa- 
tion (|lq), should behave as: 



while 



d£ P(£, t) ~ t 



-CL + P 



(20) 



(21) 
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Thus, the two above relations allow us to determine the 
exponents a and 0. The values of a„ for n = 1, ... ,6 
are shown on Fig. [l(] for p = 0.7, q = 0.1. A good fit is 
obtained for (3 = 0.48 ± 0.02 and a = 0.96 ± 0.04, in very 
good agreement with our collapsed plot. By repeating 
our analysis for other values of q (namely, 0.2, 0.3 and 
0.4), the same values for the exponents fit reasonably well 
the data. 

For q = 1/2, the different moments of the cluster 
size distribution are a\ — 0.33, a-i — 0.96, a$ — 1.60, 
a 4 = 2.22, a 5 — 2.82 and — 3.40. These exponents 
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FIG. 9. Scaling form of the cluster sizes distribution for 
p = 0.7, q = 0.1. P(£,t)t a is plotted versus £t~ 
and P = 0.5. 
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FIG. 10. Exponent a n (open circles) of the n-th moment 
of the cluster distribution function for n = 0, . . . , 6 and 
p = 0.7, q = 0.1. The line is the fit a n = 0.01 + 0.48n. 

P{£,t) ^t-^ 3 £- 4/3 , (22) 

over two decades in the variable I. Note that equa- 
tion ( p2] ) cannot be valid for arbitrary large £, because 
the moments of P(£, t) diverge with the upper limit of 
integration. 

Finally, in the domain 1/2 < q < 1, where no coars- 
ening is observed, the system approaches very rapidly its 
stationary state and no dynamical scaling has been found 
for the cluster distribution. However, in the stationary 
state, the cluster distribution takes the form: 

P(£) = Ct exp(-C 2 ^) (23) 

where C\ and C-i are two constants. 



V. INTERPRETATION OF THE RESULTS AND 
CONCLUSIONS 

The first interesting point is the particular role played 
by the value q = 1/2. As already mentioned in sec- 
tion ||, for q = 1/2, one notes the presence of an extra 
conservation law in the system. The difference between 
the average local density of particles with positive and 
negative velocities is strictly zero. It is well known that 
conservation laws has a great influence on the dynamics 
of non-equilibrium statistical systems. Accordingly, one 
may expect that the dynamics in q = 1/2 is particular. 
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FIG. 11. Scaling form of the cluster sizes distribution for 
p = 0.7, q = 0.5. P(£,t)t 1/3 is plotted versus < in a double 
logarithmic scale. The full line represents £~ 4 ' 3 . 

In view of the scaling properties of the problem it may 
be useful to think about it in terms of dynamical renor- 
malization group. Based on the results of both mean- 
field approximation and the numerical simulations, one 
is lead to conjecture the presence of three fixed points in 
this system. An unstable "critical" fixed point at q = 1/2 
and two attractive fixed points at q — and q = 1. 

When q < 1/2, the branching processes favors the ap- 
parition of pair of consecutive particles with the same 
velocities and the dynamics is governed by the attractive 
fixed point at q = 0. Large particle clusters with opposite 
velocities are formed during the time evolution and two 
kinds of interfaces are present into the system (see Fig. 

. First, let us consider the interface between two clus- 
ters of colliding particles and call this type of interface 
I 1 . Such interface has a very long life. Indeed, the prob- 
ability that a vacancy presents at one of the extremity of 
a cluster of size L traverses the cluster and perturbates 
the interface is of the order of (1 — p) L . Thus, an interface 
I\ is very stable in the long time limit where the system 
is made up of large clusters. The second type of inter- 
face, called I2, separates non-colliding clusters. Thus it 
has not necessarily a one site extension, but it can be 
wider. Hence, its behavior is more subtle. Three different 
regimes may be considered. The simplest case to discuss 
is when pa/c> 1. In this case, the interface I2 is typically 
formed by only one empty lattice site, whose dynamic is 
diffusive. Indeed, one can show that both boundaries of 
an interface I2 perform a Brownian motion. Moreover, 
when pa/c > 1, this random walk is biased, so that both 
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boundaries tend to come closer together. For sufficiently 
long time, the initial gap separating two non-colliding 
clusters will shrink to one single site, which will perform 
a random walk. Eventually, this hole will encounter an Ii 
interface, permitting the coalescence of two clusters into 
a larger one. The random walk aspect of this dynamics 
is responsible for the slow approach towards the station- 
ary state (in i -1 / 2 ) observed in the simulations. When 
pa/c = 1, the boundaries of an interface I2 both perform 
an unbiased random walk. Accordingly, the initial gap 
between two non-colliding clusters will not, on average, 
vary. However, because of the BBA dynamics, this gap 
will eventually shrink to a single site, either through a 
creation of cluster inside the gap when q 0, or through 
the coalescence of two interfaces. Thus the previous ar- 
gument holds. Finally, when pa/c < 1, the situation is 
similar: although the boundaries of I2 perform a biased 
random walk which tends to increase the separation be- 
tween the two non-colliding clusters, the coalescence of 
two interfaces or a creation of a new cluster inside the 
gap (if q ^ 0) will fill up this space in a more efficient 
way. 

Eventually, the stationary state is completely filled, 
only one cluster remains and the annihilation process do 
not act anymore. For values of q not too far from 1/2, 
this asymptotic behavior will shows up only for very long 
times. Accordingly, the results of the (finite time) numer- 
ical simulations may still be affected by the properties of 
the critical fixed point at q = 1/2, and the dynamics will 
exhibit some crossover behavior. 

In the situation q > 1/2, a majority of pairs of parti- 
cles with opposite velocities are created during branch- 
ing. Due to the annihilation processes, those particles 
will prevent the formation of large clusters of particles. 
One may anticipate that the long time dynamics is gov- 
erned by the other attractive fixed point corresponding 
to q = 1. The dynamics is no longer governed by coars- 
ening mechanism but only by the dynamics of small clus- 
ters, hence the fast (exponential like) relaxation occurs. 
Depending upon the value of p, there is a more or less im- 
portant fraction of empty sites (or holes) into the system. 
The presence of these two different dynamical regimes 
explains the jump observed in the stationary density at 
9=1/2. 

This paper shows once again, that the mean-field re- 
sults generally do not hold for low dimensional systems. 
Whereas the mean-field approximation predict the exact 
critical value for q (because the mean-field equation for ^ 
is exact when q = 1/2) and the right stationary value of 
the density when q < 1/2, it is unable to give satisfactory 
results for the density stationary value for q > 1/2, (see 
Figs [l] and |^). Unsurprisingly, the mean-field approxi- 
mation is also unable to predict the power law approach 
to the stationary state when q < 1/2, which is obviously 
governed by fluctuations. More surprisingly, its predic- 
tion of an exponentially fast approach towards the steady 
state when q > 1/2 is (qualitatively) well verified. How- 
ever, to better understand this problem, it would be use- 



ful to be able to find an exact analytical solution at least 
for the three fixed point cases (q = 0, 1/2 and 1 for arbi- 
trary values of p) as a support to the above qualitative 
picture. Unfortunately, we were not able until now to 
find such exact solutions. 

In conclusion, one sees that this simple BBA problem 
with one offspring exhibit already a very rich behavior. 
The case with two or more offsprings is a completely open 
question. 
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APPENDIX A: 

In this Appendix, we give an explicit solution to the 
mean- field equations ^ and (||) for the case q — 0. Equa- 
tions (§) and (|) lead to 



* — $ — = -c*($ 2 

dt dt v 



* 2 ). 



(Al) 



Multiplying eq. (Al) by ^ 2 , and introducing x(i) 
eq. flAA) becomes 



dx 
dt 

whose solution is 



-c*(x 2 - 1), 



x(t) 



1 + X(0) + [X(0) - 1] exp(-c Jp tt(r)rfr) 
1 + X(0) - [x(0) - 1] cxp(-c J * ¥(r)dr) 



(A2) 



(A3) 



If X(0) = ±1, thenx(i) = ±1 and one finds = ±$(t) 
for all times, where ^(i) obeys 

— =p(l T< 7*)*, (A4) 

whose solution is 

^(ty 1 =±a + exp(-pt)[*(0) _1 T <A- ( A5 ) 

Note that in these particular cases, only particles with 
velocity +c (or — c) are present in the system at all times. 
If x(0) ^ 1, one finds 



1 - 



<r$)* = x (A6) 
1 + X(0) + [x(0) - 1] cxp(-c J Q * *(r)dr 



1 + X(0) - [x(0) - 1] exp(-c JJ *(r)dr) 



(7* 







As \^(t)\ is a nondecreasing func tion of time, when depending on the sign of x(0). Thus, eq. (|A6|) reduces to 
t — » oo, the square bracket in eq. (A6) approaches ±1 eq. (A4). 



[1] Y. Elskens and H. L. Frisch, Phys. Rev. A 31, 3812 
(1985). 

[2] J. Krug and H. Spohn, Phys. Rev. A 38, 4271 (1988). 
[3] E. Ben-Nairn, S. Redner and F. Leyvraz, Phys. Rev. Lett. 

70, 1890 (1993). 
[4] P. L. Krapivsky, S. Redner and F. Leyvraz, Phys. Rev. 

E 51, 3977 (1995). 
[5] J. Piasecki, Phys. Rev. E 51, 5535 (1995). 
[6] M. Droz, P.-A. Rey, L. Frachebourg and J. Piasecki, 

Phys. Rev. Lett. 75, 160 (1995). 
[7] M. Droz, P.-A. Rey, L. Frachebourg and J. Piasecki, 

Phys. Rev. E 51, 5541 (1995). 
[8] J. Piasecki, P.-A. Rey and M. Droz, Physica A 229, 515 

(1996). 

[9] P.-A. Rey, M. Droz and J. Piasecki, Eur. J. Phys. 18, 213 



(1997) . 

[10] P.-A. Rey, M. Droz and J. Piasecki, Phys. Rev. E 57, 138 

(1998) . 

[11] J. L. Cardy and U. C. Tauber, Phys. Rev. Lett. 77, 4780 
(1996). 

[12] J. L. Cardy and U. C. Tauber, J. Stat. Phys. 90, 1 (1998). 
[13] P. Grassberger and K. Sundermeyer, Phys. Lett. 77 B, 

220 (1978); J. L. Cardy and R. L. Sugar, J. Phys. A 13, 

L423 (1980); 

[14] H. J. Kreutzer, Nonequilibrium Thermodynamics and 
its Statistical Foundations, Oxford Science Publication, 
(1981). 

[15] L. Frachebourg and P. L. Krapivsky, Phys. Rev. E 55, 
252 (1997). 



10 



